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This paper describes a method for spread- 
sheet computations of Fresnel integrals to 
six significant figures, based on successive 
improvements of known rational approxima- 
tions which are accurate to only three 
figures. Outside the range of validity of the 
improved approximations, known series 
expansions are used to obtain the Fresnel 
integrals to six figures. 
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1. Known Formulae 

The Fresnel cosine and sine integrals, 



C(x) 



dicosl— t 



spreadsheets and the algorithms were more complicated 
than necessary for this project. It was therefore decided 
to develop suitable approximations that can be 
programmed into spreadsheet programs to supply the 
required values of C(x) and S(x). 

This development was based on the following formu- 
lae listed in the NBS Handbook of Functions [3]: 



+ f(x) sinl^x 2 -g(x)cos [^ x 2 
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+ e(x) , (2a) 
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= --/(*) cos(^* 2 j -g(*)sin(^* 2 j , (lb) 
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2 + 4.142* + 3.492* 2 + 6.67* : 



+ e(x) , (2b) 



where /(*) and g(x) are auxiliary functions, are impor- 
tant in various fields of physics. This author required 
them to six significant figures for a project in optics 
which employed spreadsheet software to compute and 
plot three-dimensional diffraction patterns. It was found 
that software and algorithms for computing the Fresnel 
integrals in Fortran and C are readily available [1, 2], 
However, this software was not directly applicable in 



\e(x)\ < 2 10' 3 , 0<*<» ; (2c) 
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\e(x)\ < 3 10- 6 , |*| <0.63, (3a) 
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|e(x)|< 3 lO" 6 , |jc| < 0.73; 
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(3b) 



g(x) = g (x) + b t (x) sin[/3,(x)] 
+ b 2 (x) sin[/3 2 (x)] + b 3 (x) + s(x) 



\s(x)\ <3 \0' 6 , 0.6 < x< 3.1, 



where 



a x {x) = (1.3384 + 0.2158x)10' 3 , 



(5b) 
(5c) 

(5d) 



|e(x)|< 3 10' 6 , x>2.9 , 



(4a) 



a 2 (x) = (6.158 + 0.073x) 10~ 



a 3 (x) = (-15.604 + 2.88x)10" 



(5e) 



(5f) 



g(*) = 



15 945 



+ s(x) , 



|e(x)|< 3 10' 6 , x>2.: 



(4b) 



As the rational approximations of Eqs. (2a) and (2b) 
yield C(x) and S(x) only to ± 2.5 X 10~ 3 , it was 
necessry to improve them over a sufficiently wide range 
of intermediate arguments so that they could be supple- 
mented with the series expansions of Eqs. (3a) to (4b) to 
compute C(x) and S(x) to ±3 X 10~ 6 for arbitrary 
values of x. 



a,(x) = 
(3.0824 + 14.8294x + 0.00462x 2 )/(l+x) , (5g) 

a 2 (x) = 



(0.857 + 11.78x + 0.773x 2 )/Vl + 2.354x 2 



and 



fo,(x) = (1.367-0.2278x)10- 3 , 



(5h) 
(5i) 



2. Improvement of Rational Approxima- 
tions 

The improvement of the rational approximations of 
Eqs. (2a) and (2b) was based on a comparison of the 
numerical values of fo(x) and go(x) to the exact values 
tabulated in Ref. [1]. It was found that the residuals e(x) 
were skewed oscillatory functions that could be approx- 
imated by expressions of the form e(x) ~ A(x) = a(x) 
sin[a(x)] + b(x), where a(x) and b(x) are linear func- 
tions of x and a(x) is chosen so that the oscillations of 
e(x) are imitated. By a judicious choice of parameters it 
was arranged that the second residuals, e(x)-A(x), 
also exhibited a sinusoidal behvior so that the same 
procedure could be repeated. As the techniques used 
were empirical, based on trial and error rather than 
mathematical rigor, no further details will be given. The 
final results are: 



/CO =fo(x) + fli(x) sin[a,(x)] 
+ a 2 (x) sin[a 2 0O] + a 3 (x) + e(x) , 



(5a) 



b 2 (x) = (8.94-0. 19x)10~ 5 , (5j) 



b 3 (x) = (-7.86+11.2x)10- 5 , (5k) 



j3 t (x) = (16.442 - 0.00012x + 0.00034x 2 ) 



X (1 - e -°- 7738 ^) , (51) 



p 2 (x) = (1.524 + 29. lx - 0.392x 2 ) 



/i/1 + 5.934x 2 . (5m) 



The digits appearing in these expressions are significant 
within the uncertainties stated in Eq. (5c). By substitut- 
ing the numerical values of/(x) and g(x) given by these 
equations into Eqs. (la) and (lb), and comparing the 
results to the exact values in Ref. [3], it was found that 
they yield C(x) to ± 3.5 X 10~ 6 for 0.6 < x < 3.1. 
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3. Application 

The results of Sec. 2 were entered into a spreadsheet 
program and C(x) and S(x) were computed for 
0.6 < X < 3 from Eqs. (la) and (lb) and Eqs. (5a) to 
(5m). For the respective ranges x < 0.3, and x > 3, the 
the Taylor series of Eqs. (3a) and (3b) and the asymp- 
totic expansions of Eqs. (4a) and (4b) were used. 
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